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Abstract 

We construct explicit solutions to the discrete motion of discrete plane curves that has 
been introduced by one of the authors recently. Explicit formulas in terms the r function 
are presented. Transformation theory of the motions of both smooth and discrete curves is 
developed simultaneously. 

2010 Mathematics Subject Classification: 53A04, 37K25, 37K10, 35Q53. 

Keywords and Phrases: 

discrete curves; discrete motion; discrete potential mKdV equation; discrete integrable sys- 
tems; r function; Backlund transformation. 



1 Introduction 



Differential geometry has a close relationship with the theory of integrable systems. In fact, many 
integrable differential or difference equations arise as compatibility conditions of some geometric 
objects. For instance, it is well known that the compatibility condition of pseudospherical sur- 
faces gives rise to the sine-Gordon equation under the Chebyshev net parametrization. For more 
information on such connections we refer to a monograph [37] by Rogers and Schief. 

The above connection between the differential geometry of surfaces and integrable systems 
has been known since the nineteenth century (although the theory of integrable systems was not 
yet established). However, it is curious that the link between the differential geometry of curves 
and integrable systems has been noticed rather recently. Actually Lamb [28] and Goldstein and 
Petrich [ ] discovered an interesting connection between integrable systems and the differential 
geometry of plane curves. Namely, they found that the modified Korteweg-de Vries equation 
(mKdV equation) appears as the compatibility condition of a certain motion of plane curves. Here 
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a motion of curves means an isoperimetric time evolution of arc-length parametrized plane curves. 
More precisely, the compatibility condition implies that the curvature function of a motion should 
satisfy the mKdV equation. As a result, the angle function of the motion satisfies the potential 
modified Korteweg-de Vries equation (potential mKdV equation). 

On the other hand, in the theory of integrable systems, much attention has been paid to dis- 
cretization of integrable differential equations preserving integrability, after the pioneering work 
of Ablowitz and Ladik [ ] and Hirota [16-20]. Later, Date, Jimbo and Miwa developed a unified 
algebraic approach from the view of so-called the KP theory [5-9,27, 32]. For other approaches 
to discrete integrable systems, see, for example, [33, 38]. Thus one can expect the existence of 
discretized differential geometric objects governed by discrete integrable systems. This idea has 
been realized by the works of Bobenko and Pinkall [ ] and Doliwa [ ] where the discrete ana- 
logue of classical surface theory has been proposed, and it is now actively studied under the name 
of discrete differential geometry [4]. 

In the case of the discrete analogue of curves, Doliwa and Santini studied continuous motion 
of discrete curves in the 3 -sphere, and obtained the Ablowitz-Ladik hierarchy in [ ], where a 
semi-discrete (discrete in space variable and continuous in time variable) analogue of the mKdV 
equation is derived as the simplest case. Their formulation includes the plane curves as a limiting 
case. Hisakado et al proposed a discretization of arc-length parametrized plane curves [ ], and 
obtained another semi-discretization of the mKdV equation. Hoffmann and Kutz [25, 26] consid- 
ered discretization of the curvature function. By using their discrete curvature function and Mobius 
geometry, they obtained a semi-discrete mKdV equation that is the same as the one in [1 1]. 

The discretization of the time variable of discrete curve motion in the 3-sphere was studied 
in [12, 13], where the evolution corresponded to the discrete sine-Gordon equation [18]. As is well 
known, binomial motion of space curves induces the nonlinear Schrodinger equation via Hasimoto 
transformation [15,28]. Discretization of this curve motion has been formulated in [24,36]. 

Recently one of the authors of the present paper formulated a full discretization of the motion 
of plane discrete curves [ ] in a purely Euclidean geometric manner, where the discrete potential 
mKdV equation proposed by Hirota [ ] is deduced as the compatibility condition. It admits a 
natural continuous limit to the potential mKdV equation describing continuous motion of a smooth 
plane curves. 

In the smooth curve theory, the potential function coincides with the angle function of a curve, 
a primitive function of the curvature. However, in the discrete case, the potential function and 
the angle function become different objects. In this framework, the primal geometric object is 
the potential function rather than the curvature (see [12, 31] and Section 2 of the present paper). 
Natural and systematic construction of the discrete motion of the curves is expected by using the 
theory of discrete integrable systems. 

The purpose of this paper is to construct explicit solutions to discrete motion of discrete plane 
curves by using the theory of r functions. This paper is organized as follows. In Section 2, we 
prepare fundamental ingredients of plane curve geometry and motions (isoperimetric time evolu- 
tions) of plane curves described by the potential mKdV equation. Next we give a brief review of 
the discrete motion of discrete curves [31]. In Section 3, we shall give a construction of motions 
for both smooth and discrete curves by the theory of r functions. More precisely we introduce a 
system of bilinear equations of Hirota type, which can be obtained by a certain reduction of the 
discrete two-dimensional Toda lattice hierarchy [27,39,40]. We shall give a representation formula 
for curve motions in terms of the r function. 
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One of the central topics in classical differential geometry is the transformation theory of curves 
and surfaces. The best known example might be the Backlund transformations of pseudospherical 
surfaces. The original Backlund transformation was defined as a tangential line congruence sat- 
isfying the constant distance property and constant normal angle property (see [ ] ). In plane 
curve geometry, Backlund transformations on arc-length parametrized plane curves can be defined 
as arc-length preserving transformations satisfying the constant distance property. Such transfor- 
mations can be extended to transformations on smooth curve motions via the transformation of 
solutions to the potential mKdV equation. Motivated by this fact, we shall introduce Backlund 
transformations for discrete motion of discrete curves in Section 4 (compare with [25]). In particu- 
lar we shall give another type of Backlund transformation on motions of both smooth and discrete 
curves, which is related to the discrete sine-Gordon equation. In Section 5, we shall construct and 
exhibit some explicit solutions of curve motions, namely, the multi-soliton and multi-breather so- 
lutions. We also present some pictures of discrete motions of discrete curves. We finally give some 
explicit formulas for the Backlund transformations of both smooth and discrete curve motions via 
the t functions. 



2 Motion of plane curves 



Let y(x) be an arc-length parametrized curve in the Euclidean plane R 2 . Then the Frenet equation 
of y is 

' -K 
K 



r = 



r 



(2.1) 



Here ' denotes the differentiation with respect to x, and the function k is the curvature of y. Let us 
consider the following motion in time t, i.e., isoperimetric time evolution: 



dt 7 = 







K" + K - 

2 



-k" 

2 



7 



(2.2) 



Then the potential function 8{x, t) defined by k = 6' satisfies the potential mKdV equation [14, 28] : 

1 



o. 



(2.3) 



The function 6 is called the angle function of y in differential geometry. Note that y' can be 
expressed as 

cos 6 
sin# 



y 



For any non-zero constant A, the set of equations 



d_ 

dx 

d_ 

dt 



e + e\ 

2 

(9 + 6} 



e-e 



2A sin ■ 



6 — 6 6 — 6 

-A {(6 X ) 2 + SA 2 } sin — — + 2AB XX cos — — + 4A 2 6 X , 



(2.4) 



(2.5) 
(2.6) 
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defines a solution 6 to the potential mKdV equation [41]. The solution 9 is called a Bdcklund 
transform of 6. 



Definition 2.1 A map y 



I 2 ; n i — ► y n is said to be a discrete curve of segment length a n if 

y»+i - 7n 



Cl rl 



1. 



(2.7) 



We introduce the angle function W„ of a discrete curve y by 



A discrete curve y satisfies 



a n 

J n+l ~ 7n 



COS 

sin *F„ 



R(K n ) 



In ~ Jn-l 



a n i 
for K n = - *F n _i, where denotes the rotation matrix given by 

R(K„)-[™ K - - si "M. 
\^ sin^ n cosK n J 

Now let us recall the following discrete motion of discrete curve y™ : Z 2 
Matsuura [31]: 



y m _ y m 
I n+l in 



m _ m 
<_ n+l I n 

a n 



= 1, 



m _ m 
In I n-\ 

&n-\ 



-V m +1 _ -y m y m — y m 
in in D/'Wn "n+l "n 
7 = K k W n ) . 

o m a n 



(2.8) 
(2.9) 

(2.10) 
L 2 introduced by 

(2.11) 
(2.12) 
(2.13) 



where a n and b m are arbitrary functions of n and m, respectively. Compatibility of the system 
(2. 1 1)— (2. 13) implies the existence of the potential function &" defined by 



T i/m n ^n+l v m w n+l 

W, = j ' K " = 2 ' 

and it follows that 0™ satisfies the discrete potential mKdV equation [21]: 



tan " Urstl^tanP ^ 

\ 4 ) b m -a n \ 4 



Note that the angle function can be expressed as 



VJ/n 



m + ©'« 



(2.14) 



(2.15) 



(2.16) 



Remark 2.2 The potential discrete mKdV equation (2.15) has also been known as the superposi- 
tion formula for the modified KdV equation (2.3) [41] and the sine-Gordon equation [2,37]. 
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3 The r function representation of plane curves 

In this section, we give a representation formula for curve motions in terms of t functions. 

Let = t%(x, t; y) be a complex- valued function dependent on two discrete variables m and n, 
and three continuous variables x, t and y, which satisfies the following system of bilinear equations: 



1 9 

-D D x m -T m = - ( T * m ) 



x r„ -t 



0, 



(z>2 + a) < • t*: = o, 



_*m *m 



Dy " — a n T n +i T n ' 

^ ( n ( n _ ''w* n+1* n' 

, *m+l m _ *m m+1 / _ , x *m+l m _ q 



(3-D 
(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 



Here, * denotes the complex conjugate, and D x , D y and D, are Hirota's bilinear differential opera- 
tors (Z)-operators) defined by 

Vjyty f-g= (d x - dj{d y - d y )\d t - d f ) k f(x,y, t)g(x',y', t')\ x=x , y=yt=f . (3.7) 

For the calculus of the D-operators, we refer to [22]. In general, the functions satisfying the bilinear 
equations of Hirota type are called r functions. 

Theorem 3.1 Let t™ be a solution to equations (3.1)-(3.6). Define a real function ®™(x, t;y) and 
an R 2 -valued function y™(x, t; y) by 



&:(x,t;y) 
7':(x,t;y): 



2 r m 



1 / 

. 2V r Tr 8 ^T/v 



(3.8) 



(3.9) 



(1) For any m, n e Z and y e R, the functions 6(x, t) = ®%(x, t; y) and y(x, t) = y£(x, t\ y) satisfy 
equations (2.1)-(2.3). 

(2) For any x,t,y e R, the functions 0™ = ®™(x, t\ y) and y™ = y%(x, t; y) satisfy equations 
(2.11) (2.15). 

Proof. (1) Express y'" = '(X™,Y™). Then by using equation (3.1) together with its complex 
conjugate, we have 



PC)' = —logco^, = -\ 



In n —*m _*m In n —m _m 

2 LJ X LJy I n I n 2 LJ X LJy l n l n 



(r *m)2 



+ 



cos 0" 
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Similarly we obtain (7™)' = sin©" 1 . Differentiating (y™)' = '(cos sin®" 1 ) by x and noticing 
that k = 0' , we obtain equation (2.1): 



air = (®:y 



A - sin 0" 



-k 

cos 0™ / " \ k 
On the other hand, differentiating (y™)' by t, we have 



(TnX = (®:\ ( ^ ) = ( & n\ [l Oj 

By using the bilinear equations (3.2) and (3.3), (&") t can be rewritten as 



<-> n -r m -r* m 

Z u t ~n ~ n 

r—r T m T *m 

V I n n 

;| : HI 



o r)3 _m , x *m 
z ^ x ' n ' n 



—m T *m 



V=I 

2 



log-^ 



.* m 

n / jawe 



*m 
nix 



- Kxx r 



which yields equation (2.2). Here we have used the relation 



r\2 _m T *m , _„ 

= = (log T>*X + log I , 



1 n n 



(3.10) 



which is a consequence of equation (3.2). The potential mKdV equation (2.3) follows immediately 
from equation (3.10) by noticing that /<- = ©'. 

(2) From equation (3.4) and its complex conjugate we have 



log 



n+1 



n I y 



m m 
1 n+1 ' n 
~m _m ; 
'n+1 'n 



log 



n+1 



^n 



'n+1 'n 
1 n+1 1 n 



Adding these two equations we obtain 



(3.11) 



(logCiTT+i) v -(logr;;VT) =-a„ 



n+1 n 
'n+1 'n 



+ 



—m _jn 
'n+1 'n 
T *m — *ni 
' n+1 n 



which yields 



" +1 " =cos^', ^ = log J*tiZL = U " +1 + U » , 

" V^T Ut+l^n"/ 2 

Subtracting the second equation from the first equation in (3.1 1) we have 



ym ym 

1 n+1 L n 



sinY! 



(3.12) 



(3.13) 



Therefore we obtain 



7C +1 -y? /cos¥™ 
sin^ 



(3.14) 



which gives equation (2.1 1). Next, from equation (3.14) we see that 



y" - y'" 

I n+l / n 



y m _ m 
In I n-l 



VT»m _ mm 
1 x n -l - 



0"> _ (hV« 



(3.15) 



which is nothing but equation (2.12). Similarly, starting from equation (3.5) and its complex 
conjugate we obtain 



In In 



cosO™ 
sin 0" ! 



O m = 



1 



V=i 



log 



_m+l_m 
' n n 

-[■*m+l-r* m 

n n 



V -' It 1 11 



(3.16) 



which yields 



m+l _ m <y m _ y m 

I n fn_ _ i>((f\m _ \rjm\ I n+l in 

b ~ { n n) a 

u m u. n 



_ vpm _ 







.m+l 



0'" 

7! + l 



n n 



= W" 



(3.17) 



This is equivalent to equation (2.13). 

Finally let us derive the discrete potential mKdV equation (2.15). Dividing equation (3.6) and 



its complex conjugate by t*™ +1 t*™ +1 we have 



b m exp 



b m exp 



^®:+A 



V-L0 



.m+l \ 



a n exp 



- a n exp 



( V^e 



m+l \ 



— ifln b ln ) 



(fin ^m) 



_*m+ 1 _m 
' n+l l n 
—*m+l T *m 
1 n 1 „+l 



_m+ 1 _* m 
'«+! ' » 

_*m+l 7 -*m 
' n ' n+l 



(3.18) 



respectively. Dividing these two equations we obtain 



b m exp y-^pA - a„ exp j 

I V-f 0" ,+1 \ / V 3 ! , \ 

6 m exp 1^ 1 2 u " j - a„ exp ^ — ^ 1 



exp 



(3.19) 



which is easily verified to be equivalent to equation (2.15). Thus we have completed the proof of 
Theorem 3.1. □ 

Corollary 3.2 (Representation formula) 

The function y% can be expressed in terms of the potential function 0™ as follows: 



J cos&;Xx',t;y)dx' 
J smQ™(x', t;y)dx' 



n~l 



&Z(x,t;y) + &2 +1 (x,t;y) 




(3.20) 



Proof. The first equation is a consequence of 



d 



cos ® m (x, t;y) 
sin0™(*, t;y) 



(3.21) 
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and the second equation follows from equation (3.14). 



□ 



It should be noted here that the bilinear equations (3. 1)— (3.6) are derived from the reduction of 



the equations 

l -D x D y -C(s) ■ t» = -r^s + \K(s - 1), (3.22) 

(D 2 X -D Z ) ^{s + \)--C(s) = Q, (3.23) 

(d 3 x + D t + 3D X D Z ) r';;(s + 1) • t™(s) = 0, (3.24) 

D y CiW ' <(s) = -a n ^ +1 (s + \K(s - 1), (3.25) 

D y ^ n n + \s) ■ r» = -b m -C +l (s + \K(s - 1), (3.26) 

b m -C\s + DCiW - amt x {s + 1K + V) + (a n - b m yCti? + = 0, (3.27) 

for t™(s) = t™{x,z, t;y; s), which are included in the discrete two-dimensional Toda lattice hierar- 
chy [27, 39,40]. In fact, imposing the conditions 

--C(s) = £<(*), t';:( S + 1) = Ct*:( S ), B,CeR, (3.28) 



and denoting t% = r™(0), then equations (3.22)-(3.27) yield equations (3. 1)— (3.6), respectively. 



4 Backlund transformations 



We start with the following fundamental fact on plane curves. 

Proposition 4.1 Let y(x) be an arc-length parametrized curve with angle function 9{x). Take a 
non-zero constant A and a solution 0(x) to 



(e + e\ 



= 2A sin ■ 



e-e 



(4.1) 



Then 



y(x) = y(x) + jR 



( 9{x) - 6{x) 



y'(x) 



(4.2) 



is an arc-length parametrized curve with angle function 9(x). In other words, ify(x) is a solution 
to equation (2.1), theny(x) is another solution to equation (2.1) withH(x) = 6'(x). The curve y is 
called a Backlund transform ofy. 

Proposition 4. 1 can be verified easily by direct computation. We next extend the Backlund trans- 
formation to those of the motion of a curve. 

Proposition 4.2 Let y(x, t) be a motion of an arc-length parametrized curve determined by equa- 
tions (2.2) and (2.3). Take a Backlund transform 6(x,t) defined by equations (2.5) and (2.6) of 
6(x, t). Then 



1 

y(x, t) = y(x, t) + -R 
A 



(6(x, t)-9(x, t)\ 



y'(x,t) 



(4.3) 



is a motion of an arc-length parametrized curve with the angle function 9(x, t). 



8 



Proof. By the preceding proposition, y satisfies the isoperimetric condition \y'\ = 1. Computing 
the /-derivative of y by using (2.6), we can show that y satisfies equation (2.2) with 7c = 6' □ 

Now we introduce a Backlund transformation of a discrete curve. 

Proposition 4.3 Let y n be a discrete curve of segment length a n . Let 0„ be the potential function 
defined by 



Jn+l ~ Jn 



a„ 



cos *F„ 
sin 



For a non-zero constant A, take a solution 0„ to the following equation: 



tan 



©n+l - ®« 



tan 



®n ~ ®n 



+ 1 



Then 



1 

7n = Jn + /? 



f0„-©„ + l^ 



7n+l ~ Jn 



a n 



is a discrete curve with the potential function 0„. 
Proof. It suffices to show that 

7n+l -Jn 



a„ 



cos W n 
sin 



~ _ 0„ + i + e n 



for y n defined by equation (4.6). This follows from equations (4.4) and (4.5). □ 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



We next extend the Backlund transformation to those of the motion of a discrete curve. In order 
to do so, we first present the Backlund transformation to the discrete potential mKdV equation. 

Lemma 4.4 Let 0™ be a solution to the discrete potential mKdV equation (2.15). A function &" 
satisfying the following system of equations 



tan 



tan 



( Q m +1 







m \ 



tan 



( Gfi _ f\m \ 



tan 



( (5)'" _ (5) 

a 



m+l \ 



n 



(4.8) 
(4.9) 



gives another solution to equation (2.15). We call &" a Backlund transform of®™. 
Proof. First note that equation (2.15) is equivalent to 



U m+l +U" 

g n T w n 



Tjm+l.Tjm a n I jjm 

e u n+i * u n+i = — [e "+ 

h V 



IJIH+ 1 ±Jjm+ 1 

e ■ "+ l : 



■m+ \ \ 



(4.10) 



V— 10'" 

where we put — ^ = U™ for notational simplicity. Similarly, equations (4.8) and (4.9) are 
rewritten as 



jim ,jjm TJ m+ ' 4-TI m+ ^ 

g n n g n n 



e u ^ u ^=Aa n (e u ^ +u "-e u ^ +u "), 
Ab m (e u « +l+u » -e U » +l+U »), 



(4.11) 
(4.12) 
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respectively, where — j -2 " = U%. Subtracting equation (4.12) from equation (4.11), we have 

e ur+ur _ e u: +1 +K +1 = A ^ ne ^ +1 *K _ hm e^^) _ A ^ a j:^;; _ bme u '" +> ^ (4 13) 

Similarly, subtracting equation (4.12)„^ n+1 from equation (4.1 l) m _ >m+1 , we get 

fjm+l.jjm+\ Tim .Tim ( ijm+l.rim+l jjm+\,jjm \ / Tjm+l + Tjm+l Tjm+1 + fjm \ 

e u n +u„ _ e u„ +l +u n+l _ x\a n e u »^ +u " - b m e u ^ +ty "+iJ - A [a n e u »^ +u " -b m e u ^ +u ^). (4.14) 

Subtracting equation (4.14) from equation (4.13) yields 

/ J]'" jSj m u" ,+l +u m+, \ i I u m+l +ij m u m+l +u'" \ 

v ' v ' (4 15) 

/ J]"' + f/'» U"' + ' + U"' +1 \ i ( u m+ ' + U'" U"' + ' + U"' \ 

= a n ie ■<+< "- e ><+> " \-b m [e " +u >< - e »+> ><+i\ . 

Now we see that the right-hand side of equation (4.15) vanishes since it is equivalent to equation 
(4. 10). Then the left-hand side gives equation (2. 15) for 0™. □. 



Proposition 4.5 Let y™ be a discrete motion of a discrete curve. Take a Bdcklund transform O™ of 
0J" defined in Lemma 4.4. Then 



Y n = y m + -R 

in in A 



^n ^n+l 



m _ m 
I n+l i n 

a n 



(4.16) 



is a discrete motion of a discrete curve with potential function 0™. We call y™ a Bdcklund transform 

ofrt- 

Proof. It suffices to show that satisfies equations (2.1 1)— (2. 13) with potential function 0™. But 
equations (2.1 1) and (2.12) follow from Proposition 4.3 immediately. Noticing the symmetry in n 
and m, similar calculations to those in Proposition 4.3 yield 



~m+l _~m 
in in 



COS 



sin 



2 



(4.17) 



by using equation (4.9). Comparing equations (4.7) and (4.17) we obtain 



~m+l _~m 
in in 



R 



( 0'«+i 







~m _~m 
' n+l I n 

a n 



(4.18) 



which implies equation (2.13). □ 



It is possible to construct another type of Backlund transformation for motions of both smooth 
and discrete curves by using the symmetry of the potential mKdV equation (2.3) and discrete 
potential mKdV equation (2.15). In fact, if 8(x, t) is a solution to equation (2.3), then -8{x, t) 
satisfies the same equation. Combining this symmetry and the Backlund transformation defined 
by equations (2.5) and (2.6), we have the following Backlund transformation. 
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Lemma 4.6 Let 6{x,t) be a solution to the potential mKdV equation (2.3). For any non-zero 
constant A, a function 6(x, t) satisfying the following set of equations 



d (e-e\ 
d_ (e-e\ 

dt 



2A sin ■ 



e + e 



= -A {(6 X ) 2 + SA 2 ) sin ^— - 2A6 XX cos ^ - AA 2 9 X , 



e + e 



(4.19) 
(4.20) 



gives another solution to equation (2.3). 

Lemma 4.6 immediately yields the following Backlund transformation for y(x, t). 

Proposition 4.7 Let y(x, t) be a motion of an arc-length parametrized curve determined by equa- 
tions (2.2) and (2.3). Take a Backlund transform e(x, t) of e(x, f) defined in Lemma 4.6. Then 



y(x, t) = S 



y(x, t) + jR 



e(x, t) + e(x, t) 



Y(x, t) 



s = 



1 
-1 



(4.21) 



is a motion of an arc-length parametrized curve with angle function e(x, t). 

Note that equations (4.19) and (4.20) can be derived from equations (2.5) and (2.6) simply by 
putting e(x, t) = ~e(x, t). Moreover, putting 



y(x,t) :=y(x,t) + -R 
A 



e(x, t) + e(x, t) 



y'(x, t), 



and noticing equation (2.4) and Proposition 4.2, we have 

y(x, t) = 



cos(~e(x, t)) 
sin(-#(;t, 0) 



cos e(x, t) 
- sin #(.*:, t) 



(4.22) 



which implies Proposition 4.7. 

Similarly, if 0™ is a solution to equation (2.15), then -0™ satisfies the same equation. Therefore 
Lemma 4.4 and Proposition 4.5 lead to the following Backlund transformations. 

Lemma 4.8 Let 0J" be a solution to the discrete potential mKdV equation (2.15). A function 0™ 
satisfying the following system of equations 



tan 



tan 



( c-)'" + C-Y" \ 



tan 



A a " 



( ©m+l + @» 



\ + b m 

~~ - - h 



tan 



( + (5)^+1 \ 



(4.23) 
(4.24) 



gives another solution to equation (2.15). 
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Proposition 4.9 Let y™ be a discrete motion of a discrete curve. Take a Bdcklund transform 0™ of 
0™ defined in Lemma 4.8. Then 



y m _ m 

a n 



S = 



1 
-1 



(4.25) 



15 a discrete motion of a discrete curve with potential function 0™. 
Remark 4.10 

(1) It may be interesting to point out that equation (4.23) and equation (4.24) can be rewritten as 



sin 



sin 



w «+i w « ^ w r 



/la„ sin 



= Ab m sin 



v>> n 



(4.26) 
(4.27) 



respectively, which are essentially equivalent to the discrete sine-Gordon equation [18]. 

(2) The Backhand transformations described in Propositions 4.2 and 4.5 satisfy the 'constant 
distance property', i.e., |y- y\ = l/A or |y™ - y'"\ = l/A. These transformations may be 
regarded as the one-dimensional analogue of the original Backlund transformations of the 
pseudospherical surface [37]. On the other hand, the Backlund transformations proposed in 
Propositions 4.7 and 4.9 are characterized by the property \y-Sy\ = l/A. 



5 Explicit solutions 
5.1 Solitons and breathers 

For N e Z> we define a function t%(s) = t%(x, t; y, z; s) by 



r^(5) = exp 



/ n-l m-1 \ 



V n' 



m' ' 



for (x, t;y, z) e R 4 and (m,n,s) e Z 3 . Here / s w = f s ,(x, t; y, z', m, n) (i = 1,...,N) satisfies the 
following linear equations: 



=1,...JV 



(5.1) 



Of 



At) 



fi\m, n) - jf(m, n - 1) (0 /j !) (m, «) - /^(m - 1, n) {i) 

= f s lM, n), = /; + i(m, n). (5.3) 



K0/ 



-4/; 



(0 
+3' 



dft 



(0 



r(0 



(5.2) 



o n -i 



'm-1 



12 



For N = 0,we set detC/^.j)^^ = 1. A typical example for /j W is given by 

n-l m-l 

= <*iPi - a "'P^ f[ (1 - b m 'PiT 1 e PiX+ ^- 4p ' t+ ^, 

n' m' 

n-l m-l 



(5.4) 



(5.5) 



n' 



m' 



where p h q h a t and /?, are arbitrary complex constants. 

We note that r™ and are functions of continuous variables x, y, z, t and discrete variables 
m, n, s, but we will indicate only the relevant variables according to the context, for notational 
simplicity. Then it is well known that t%(s) satisfies the bilinear equations (3.22)-(3.27) [22,27, 
29,30,34,35,40]. Actually by using the linear relations (5.2) and (5.3), equations (3.22)-(3.27) 
are reduced to the Pliicker relations, which are quadratic identities of determinants. 

It is possible to construct the solutions to the bilinear equations (3. 1)— (3.6) by imposing the 
reduction condition (3.28) on t%(s) in equation (5.1). Those conditions are realized by putting 
restrictions on the parameters of the solutions. As an example, we present the multi-soliton and 
multi-breather solutions. 



Proposition 5.1 Consider the r function 

( n-l 

t% = exp 



m-l \ 



=1,...JV 



x + 2_j ««' + ^ b n > y det (fj'\). ._ 

n' m' / _ 

fit) = e m + g« 

m-l 

*> = atf f[(l - a n , Pi T l f[d " bntPd- 1 ^"*', 

n' m' 

n-l m-l 

=M-pd' n^ 1 +a *p^ n (i + b m 'PiT i e- pix+Ap >^. 



71-1 



(5.6) 
(5.7) 

(5.8) 



(1) Choosing the parameters as 

p u a, eR, Pi G V^lR (i=l,...,N), 



(5.9) 



then T m satisfies the bilinear equations (3.1 )-(3.6). This gives the N-soliton solution to 
equations (2.3) and (2.15). 



(2) Taking N = 2M, and choosing the parameters as 

p h a h Pi G C (i = 1, . . . , 2M), p 2k = p* 2k -\ (k=l,...,M), 
azk = fok = -Pik-i ( k= 1 > • • • > M )> 



(5.10) 



then r m satisfies the bilinear equations (3.1)-(3.6). This gives the M -breather solution to 
equations (2.3) and (2.15). 
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Proof. It is sufficient to show that the conditions in equation (3.28) are satisfied. We first impose 
the two-periodicity in s, i.e., t™{s + 2) = const, x t%(s). For t%(s) in equation (5.1) with entries 
given by equations (5.4) and (5.5), putting 



we have 

which implies that 
Note that the condition 



<7, = ~Pi, (5.11) 
fS 2 = Pifl^ (5-12) 



N 



t™{s + 2)=A n t™{s), A n = Y\pI (5.13) 



dz 



i=\ 



N 



= B N ^(s), B N = Y,PI (5-14) 



is also satisfied simultaneously. Now we consider cases (1) and (2) separately. 

Case (1). We see from equations (5.7) and (5.8) together with equation (5.9) that 

jf=Pijf* (5.15) 

and so 

N 

<(l) = C«O0), C N = Y\ Pi eR. (5.16) 
Case (2). We see from equations (5.7) and (5.8) together with equation (5.10) that 

f( 2k ) — „* A2k-\)* A2k-l) _ * A2k)* ( c n , 

J\ -Pik-xjQ ' J\ - Pik Jo ' y?- l/ ) 

and so 

M 

7^(1) = C N t^(0), C n = (-1) M Y] \p 2i \ 2 6 R. (5.18) 

!=1 

Therefore we have verified that the conditions in equation (3.28) are satisfied for both cases. Then 
putting t™ = t™(0), we obtain the desired result. □ 

We present some pictures of the motions of the discrete curves. Figure 1 shows the simplest 
example of a curve, which corresponds to the 1-soliton solution (loop soliton). The next example 
illustrated in Figure 2 describes the interaction of two loops, which corresponds to the 2-soliton 
solution. Figures 3 and 4 show the motions which correspond to the 1 -breather and 2-breather 
solutions, respectively. 
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Figure 1. Parameters in equations (5.6), (5.7) and (5.8): N = 1, x = 0, y = 0, a x = -l,/?i = V-T, 
P\ = 0.3, a n = l, b m = 0.5. 





m = -46 



m = -20 





m = - 1 



m = 30 



Figure 2. Parameters in equations (5.6), (5.7) and (5.8): /V = 2, x = 0, y = 0, a\ = -1, a 2 = 1, 
^ = -£ 2 = V^T, /»i = 0.3, p 2 = 0.9, a„ = 1, b m = 0.5. 

5.2 Solutions via Backlund transformations 

In the theory of integrable systems, the Backlund transformations are obtained from the shift of 
a certain discrete independent variable, which also applies to our geometric transformations. We 
first introduce discrete variables k, I, and regard the determinant size /V as an additional discrete 
variable. We then extend the r function as T%(k, I, N) = t™(x, t; y; k, I, N) in the following way: 



n' m' k' I' 1 



n-\ m-1 k-l l-l 



T™(k,l,N) = exp - x + 



'■/ i v ' 



(5.19) 



f(J) = flt +e H 3 



(5.20) 



15 



m = -1 m = 20 

Figure 3. Parameters in equations (5.6), (5.7) and (5.8): N = 2, x = 0, y = 0, a\ = a* 2 = 1, 
A = -B\ = l, Pl =p* 2 = 0.2 - 0.2 V-T, a„ = 1, b m = 1.5. 

n-l m—1 k-l l-l i 

e"' = atp\ Y\(l - a niPi y x f](l - b m , Pi y ] ]~](1 - c^pO" 1 \\ 1 " ~ 



(5.21) 



n' m' k' V 

n-l m-1 k-l l-l I , v-1 

=/3 i (-p i yY\a + an'PiT 1 Y\(i+b m , Pi y l Y\(i +c k , Pi y 1 y\ i + - 

-Pj.x+4p 3 t-—y 

x e • p i . 

Accordingly, we extend the relevant dependent variables such as and y in the same way. 
Proposition 5.2 

(1) For any k e Z, y(x, t) = y(x, t;k+ I) is a Bdcklund transform ofy(x, t) = y(x, t; k) related by 
equation (4.3) with A = —. 

c k 

(2) For any k e Z,y™ = y™(k + 1) is a Bdcklund transform ofy™ = y™(k) related by equation 
(4.16) with A = 

(3) For any N e Z> , y(x, t) = y(x, t; N + 1) is a Bdcklund transform ofy(x, t) = y(x, t; N) related 
by equation (4.3) with A = - P n+i- 

(4) For any N e Z> , y™ = y™(N +1) is a Bdcklund transform ofy™ = y™(N) related by equation 
(4.16) with A = - P n+i- 

(5) For any I e Z, y(x, t) = y(x, t\ I + 1) is a Bdcklund transform ofy(x, t) = y(x, t; I) related by 
equation (4.21) with A = dj. 

(6) For any I e Z, = y%(l + 1) is a Bdcklund transform ofy™ = y™(l) related by equation 
(4.25) with A = d h 
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m = -1 m = 20 

Figure 4. Parameters in equations (5.6), (5.7) and (5.8): N = 4, x = 0, y = 0, a\ = a* 2 = 1, a 3 = 

<= v^a = -£* = i,/3 3 = -/3* 4 = v=r,pi = ^ = 0.2-0.2^,^ = /»: = o.8 +o.8 v=r, 

^rt 1? 1.5. 

Proof. We first prove (1) and (2). It follows from equation (3.4) that the r function satisfies the 
bilinear equation 

D y r>;;(k + 1) • r:\k) = -c k T™(k + l)T™(k), (5.22) 

because of the symmetry with respect to the discrete variables m, n and k in equations (5.19)— 
(5.21). Then by an argument similar to that in the proof of Theorem 3.1, we see that 

( , / e(*+i)+g(fe) \ ^ 



y(k+\)-y(k) 

Ck 



cos 
sin 



in(M) 



(5.23) 



From equation (2.4), we have equation (4.3) with y = y(k + 1) and 6 = 6{k +1): 



**+'>-**> = fi(«* +1 >-«»W (5.24) 
Similarly from equation (3.14), we obtain equation (4.16) with = y™(k+ 1) and 0™ = ®™(k+ 1): 

y?(* + i) - )OT _ ^ / 0y(*+i)-0* +1 (*) \ r? +1 (*) -??(*) (5 25) 



Cjt \ 2 / a„ 

which proves (1) and (2). The statements (3)-(4) and (5)-(6) can be proved in much the same way 
as (l)-(2), by using the bilinear equations 

D y r^(N + 1) • <(A0 = —OiV)OiV + 1), (5.26) 

Pn+i 

D y T%(1 + 1) • OZ) = ~t™{1 + 1K(/), (5.27) 
respectively. These bilinear equations will be proved in Appendix. □ 



17 



Remark 5.3 Here we give a physical interpretation of the Backlund transformations described 
above. The Backlund transforms in (l)-(2) and (5)-(6) of Proposition 5.2 correspond to changing 
the phase of solitons (loops), in other words, the positions of solitons. On the other hand, the 
Backlund transforms in (3)-(4) correspond to increasing the number of solitons (loops). 

Computing the potential functions of the Backlund transforms of the curves, one can verify the 
following result. 

Corollary 5.4 

(1) For any k € Z, 6(x, t) = 6(x, t;k + 1) is a Backlund transform of9(x, t) = 6(x, t; k) related by 
equations (2.5) and (2.6) with A = — . 

(2) For any k e Z, ® m = &"(k + 1) is a Backlund transform of®™ = ® m (k) related by equations 
(4.8) and (4.9) with A = ± 

(3) For any N e Z> , 9(x, t) = 6(x, t\ N + 1) is a Backlund transform of6(x, t) = 9{x, t\ N) related 
by equations (2.5) and (2.6) with A = -p^+i- 

(4) For any N e Z> , 0^ = ®™(N + 1) is a Backlund transform of®™ = ® m (N) related by 
equations (4.8) and (4.9) with A = -p N+ \. 

(5) For any I € Z, 6(x, t) = 6(x, t;l + I) is a Backlund transform of Q(x, t) = 9(x, t\ I) related by 
equations (4.19) and (4.20) with A = d[. 

(6) For any I G Z, ® m = ®™(l + 1) is a Backlund transform of®™ = 0J"(/) related by equations 
(4.23) and (4.24) with A = d h 



A Derivation of bilinear equations (5.26) and (5.27) 

In this appendix, we show that the r function given in equations (5.19)— (5.21) actually satisfies 
the bilinear equations (5.26) and (5.27). For this purpose, we first introduce the generic r function 
r>:(k, l, N; s) = t>:(x, t; y, z; k, I, N; s) by 



T%(k,l,N; s) = exp 



r n-l m-1 k-l l-l ^ \ 1 

V n> m' k> I' Ul J J 



for (x, t; y, z) e K. 4 , (m,n,k,l,s) e Z 5 and 6 Z> . We require f^ = fjjr(x,t;y,z',m,n;k,l,N) 
(i = 1, . . . , AO to satisfy the linear equations (5.2), (5.3) and 



f ( J\k,l)-f^\k-l,l) ( , & 



C*-l 

■CO w 



d, 



i- 1 



(A.2) 



A typical example for f s is given by 



fit) = e m + ^ s 



(A.3) 
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n-\ 



m-l 



k-l 



I- 1 



n-\ 



m— 1 



k' 
k-l 



V 
l-l 



= a iP \ Y\(l - a n , Pi T l Y\(l - b m , Pi T l f[(l - cpd- 1 f[ ( 1 " j\ e p - x ~ A ^ y , 



eV = frq* Y\(l - a n , qi T l - b^qd' 1 " c ^d~ l Hi 1- - 

n' m' k' V ' 

where g,-, a, and are arbitrary complex constants. We put 

<(y;fe,/,iV;,) = det(/®._ 1 ).. =i ^. 



(A.4) 



(A.5) 



Proposition A.l The function cr satisfies the following bilinear equations: 
D y <(iV + 1; 5) • 5) = <(iV; 5 + l)<(iV + 1; s - 1), 



(A.6) 



A - <(/ + 1; s) ■ 5 + 1) = + 1; 5 + 1)<(Z; 5). (A.7) 

We apply the determinantal technique in order to prove Proposition A. 1 . The bilinear equations are 
reduced to the Pliicker relations, which are quadratic identities of determinants whose columns are 
appropriately shifted. To this end, we construct such formulas that express the determinants in the 
Pliicker relations in terms of the derivative or shift of a discrete variable of cr™(k, I, N; s) by using 
the linear relations of the entries. For the details of the technique, we refer to [22, 29, 30, 34, 35]. 
We introduce the notation 



o^(l,N,s) = \0 h l h ... , N-2 h N~h\, 
where 'j/' denotes the column vector 



Ji 



(A.8) 



(A.9) 



Lemma A.2 The following formulas hold: 

d y o^(l,N,s) = \ -1, 1, ••• , N-2, N-l\, 
o^(l+l,N,s) = \0 l+1 , 1, ■•• , N-2, N-l\, 
d l o J Z(l+\,N,s) = \\ M , 1, ••• , N-2, N-l\, 

-(d l d y -l)o J 2(l+l,N,s) = \0, 2, ••• , N-2, N-l 



(A. 10) 
(A.ll) 
(A.12) 

(A.13) 



Note that the subscripts of column vectors are shown only when I is shifted for notational simplicity. 
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Proof. Equation (A. 10) can be verified by direct calculation by using the fourth equation in (5.2). 
We have 

ct£(/+1,#;j) = |0/ + i, ••• , N-2 1+l , N-1 M \. (A. 14) 

Adding the (N- l)th column multiplied by d\ to the Nth column and using equation (A. 2), we have 

a^(l+l,N;s) = \0 M , l l+ u ■•• , A^-2 /+1 , N-l,\. (A.15) 

Similarly, adding the (i - l)th column multiplied by d\ to the zth column and using equation (A.2) 
for i = N - 1 , . . . , 2, we obtain 



a^(l+l,N;s) = \0 l+1 , 1, 



, iV-2, iV-1 |, 



(A.16) 



which is equation (A.l 1). Multiplying d\ to the first column of equation (A.l 1) and using equation 
(A.2), we obtain equation (A. 12). Finally, differentiating equation (A. 12) with respect to y yields 



d l d y u^(l+\,N;s) = \0 l+u 1, 2, ••• , N-2, N-\ \ + \ 1 M , 0, 2, 



= c^(l+\,N;s)-\0, l l+u 2, ••• , N-2, N-l |, 

which is equivalent to equation (A. 13). This completes the proof. □ 

Proof of Proposition A.l Consider the Pliicker relation (see, for example, [35]), 

0=| - 1,0,1, -- ,iV-2|x|l,--- ,N-2,N-l,<f>\ 
+ 10,1, ••• ,N-2,N-1 | x | -1,1,--- ,iV-2,0| 
-|0,1, ••• ,N-2,<f>\x\ - 1,1,--- ,N-2,N- 1 |, 



N-2, N- 1 1 

(A.17) 



(A. 18) 



where <p is a column vector given by 






1 



(A.19) 



By using equations (A. 8) and (A. 10), expanding the determinant with respect to the column 0, 
equation (A. 18) can be rewritten as 

= <(Af; s - 1) <(iV - 1; s + 1) + <(iV; s) d y a^{N -l;s)- <(iV - 1; j) <9 y <(iV; s), (A.20) 

which implies equation (A. 6). Similarly, applying Lemma A.2 on the Pliicker relation 



we obtain 



= | - 1,0, 1,- • ,N-2 | x ] /+ i, I,--- ,N-2,N-1\ 
-|0 /+ i,0,l,...,JV-2|x| -1,1,--- ,N-2,N-1\ 
— I 0, 1, - - - ,N-2,N-1 | x | -1,0 /+I ,1,-- - ,N-2\, 



= of (I; s - 1) x tj%(l +l;s)- d,o^(l + 1; 5 - 1) x 5 V <(Z; 5) 
- <(/; 5) x [-(JA, - !)<(/ + 1; 5 - 1)] , 



(A.21) 



(A.22) 
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which is equivalent to equation (A.7). This completes the proof. 



□ 



From Proposition A.l and equation (A.l), we see that T%(k, l,N; s) satisfies 



D y ^(N + 1; s) ■ t™(N; s) = ^(N; s + 1K(N + 1; s - 1), 



(A.23) 



D y t?(Z + 1; s) ■ t>;;(1; s + 1) = + 1; s + !)<(/; s). 



(A.24) 



We finally obtain equations (5.26) and (5.27) from equations (A.23) and (A.24), respectively, by 
imposing the reduction condition (3.28). 
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